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ABSTRACT 

We present XMM-Newton and Chandra X-ray observations of the middle-aged radio pulsar 
PSR B0355+54. Our X-ray observations reveal emission not only from the pulsar itself, but also 
from a compact diffuse component extending ~ 50" in the opposite direction to the pulsar's proper 
motion. There is also evidence for the presence of fainter diffuse emission extending ~ 5' from the 
point source. The compact diffuse feature is well-fitted with a power-law, the index of which is consis- 
tent with the values found for other pulsar wind nebulae. The morphology of the diffuse component 
is similar to the ram-pressure confined pulsar wind nebulae detected for other sources. The X-ray 
emission from the pulsar itself is described well by a thermal plus power-law fit, with the thermal 
emission most likely originating in a hot polar cap. 

Subject headings: pulsars: individual (PSR B0355+54) - stars: neutron - X-rays: stars 



1. INTRODUCTION 

Isolated pulsars constitute one of the most power- 
ful laboratories for studying particle acceleration in 
astrophysics. A significant fraction of the energy 
from rotation-power ed pulsars is converted into a wind 
(|R ees k, Gunnlll974|) . which travels at a velocity close to 
the speed of light. The interaction of this pulsar wind 
with the ambient medium produces a shock and acceler- 
ation of the relativistic particles at the shock generates 
synchrotron emission. This non-thermal diffuse emission 
manifests itself as a pulsar wind ne bulae (PWN) or vle- 
rion at radio an d X-ray energies (e.g. lRees fc Gunnll974t 
lGaensleill2001|) . Due to the short synchrotron lifetimes 
of high energy electrons, X-ray emission from a PWN 
directly traces the current energetics of the pulsar. The 
spectral and morphological characteristics of an X-ray 
PWN therefore reveal the structure and composition of 
the pulsar wind and the orientation of the pulsar's spin 
axis and/or velocity vector. 

The middle-aged 156 ms radio pulsar PSR B0355+54 
is known to emit X -rays ffielfandlll98llSeward fc Wand 
Il988t ISland Il99l and gamma-ravs l|Bhat et alJ Il990j) . 
Helfand (1983) reported the first detection in X-rays 
of the source using data from Einstein, stating that 
emission extended 5' from the pulsar. However, 
ISeward &: Wand l|1988(l analyzed the Einstein data and 
concluded that while there was evidence for weak emis- 
sion 1.7' from the source position, emission from the pul- 
sar itself was not detected. Nevertheless, they did not 
rule out the possib ility that the emission could be associ- 
ated with a PWN. Iila^lll99l detected PSR B0355+54 
in a 20 ks ROSAT observation, but owing to the lack of 
counts it was not feasible to perform a spectral analysis. 
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The analysis of the ROSAT data also led to the detection 
of faint e xtended emi ssion ~ 1.6' from the pulsar posi- 
tion, but Slanej l|1994|) did not believe there was enough 
evidence to support a link between the source and the 
extended emission. 

In this paper, we report on XMM-Newton and Chan- 
dra observations of PSR B0355+54 which we use to in- 
vestigate the presence of diffuse emission that can be 
attributed to a PWN. 

2. OBSERVATIONS AND DATA REDUCTION 

PSR B0355+54 was observed with XMM-Newton on 
2002 February 10 for 29 ks. We used data from the 
Europ ean Photon Imaging Camera (EPIC) PN instru- 
ment (Striidcr et alJl2001l) for the spatial, spectral and 
timing analysis. The PN was configured in small win- 
dow mode and the thin blockin g filter was us ed. Data 
from the MOS1 instrument ijTurner et al.ll200lTl was also 
used for the spatial analysis. The MOS1 was operated 
in full window mode with the medium filter. We reduced 
the EPIC data with the XMM-Newton Science Analysis 
System (SAS version 6.1.0). In order to maximize the 
signal-to-noise ratio for our XMM-Newton observation, 
wc filtered the data to include only single, double, triple 
and quadruple photon events for the MOS1, and only 
single and double photon events for the PN. Data were 
filtered to exclude events that may be incorrect, for ex- 
ample those next to the edges of the CCDs and next to 
bad pixels. We only included photons with energies in 
the range 0.3 — 10 keV. 

PSR B0355+54 was observed for 66 ks on 2004 July 
16 with the ACIS-S array on Chandra in the very faint 
timed exposure imaging mode. We performed standard 
data processing using CIAO version 3.2. The data were 
filtered to restrict the energy range to 0.3 — 10 keV and 
to exclude times of high background. 

3. SPATIAL ANALYSIS 

Initial inspection of the images created from the EPIC- 
PN and EPIC-MOS1 data show relatively strong emis- 
sion at the pulsar position and evidence for extended 
emission near to PSR B0355+54 (see Figure ^ top 
panel). We generated a mosaic of the PN and MOS1 im- 
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Fig. 1.— X-ray detection of PSR B0355+54 and its diffuse 
emission. Top panel: Gray-scale plot of the 0.3 — 10 keV XMM- 
Newton PN image. Middle panel: Gray-scale plot of the 0.3 — 10 
keV Chandra AGlS image. The boxes define the regions used to 
extract spectra for the diffuse emission. Bottom panel: The same 
image as the middle panel with the contribution from the X-ray 
point source removed. The image is smoothed with a Gaussian of 
width ~ 2". The arrow shows the direction of the proper motion 
of the pulsar and has a length of 20" . The circle marks the "south 
west" source (see Section 3). 



ages and measured the X-ray source positions using the 
SAS source detection tool EDETECT.CHAIN. We com- 
pared the positions of the field stars in our observation 
with the positions from optical catalogs to determine an 
astrometric correction. This correction was applied to 
the X-ray coordinates of the pulsar, resulting in R. A. 
= 03 h 58 m 53?73, decl. = +54°13'12'/12 (J2000), with an 
rms error of Of! 78. This position lies l'/6 from the radio 
position. 

To confirm the presence and examine the extent of the 
diffuse emission in the XMM-Newton data we have com- 
pared the detected PN emission with that for a point 
source. We calculated the intensity for the pulsar by 
using bilinear interpolation at regularly spaced points 
along the direction of p roper motion of PSR B0355+54 
ijChatteriee et al.ll2004|) . We compared this profile with 
the XMM-Newton point-spread function (PSF) for the 
PN at 1.5 keV, which we generated using the King profile 
parameters included in the XMM-Newton calibration file 
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Fig. 2. — X-ray emission from PSR B0355+54 as a function of 
distance from the point source along the direction of the proper 
motion of the pulsar (solid line), compared to the instrument PSF 
determined at 1.5 keV and the location of the pulsar (dashed line). 
Top panel: XMM-Newton PN. Bottom panel: Chandra ACIS. 



"XRT3_XPSF_0006.CCF.plt" 6 . In Figure H (top panel) 
we show the profiles for the pulsar and the PN PSF. 

The positions of the X-ray sources for the ACIS data 
were determined using the CIAO source detection tool 
WAVDETECT. The image does not contain enough 
sources with known counterparts to perform an astro- 
metric correction to the coordinates. A point source is 
detected at R. A. = 03 h 58 m 53f70, decl. = +54°13T3'.'87 
(J2000), which is 0'.'14 away from the radio pulsar po- 
sition. This source is consistent with being the X-ray 
counterpart of the pulsar. 

The ACIS image also reveals a faint tail of emission 
in the opposite direction to the pulsar's proper motion 
(see Figure n middle and bottom panels). Again wc 
determined the net counts from the source and diffuse 
emission at regularly spaced intervals along the direc- 
tion of the pulsar's proper motion. We generated a 
PSF for PSR B0355+54 using the Chandra PSF library 
evaluated at 1.5 keV and the location relevant to our 
source. The PSF was normalized to the total counts in 
PSR B0355+54. We calculated the net counts for the 
PSF in the same intervals as for PSR B0355+54. The 
source and PSF profiles are shown in Figure |3 (bottom 

6 See http: / / xmm.vilspa.esa.es / docs/documents/CAL-SRN-0100-0-0.ps.gz 
for more information 
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Fig. 3. — XMM-Newton MOS1 X-ray emission from 
PSR B0355+54 as a function of distance from the point source, 
orientated 20° further North from the direction of the pulsar's 
proper motion (solid line), compared to the mean background (dot- 
ted line). 



panel). 

The XMM-Newton PN image shows emission ~ 45" 
south west of the pulsar which could be a source, how- 
ever, there is no corresponding emission at this position 
in the Chandra ACIS image (see Figure ^ . We per- 
formed wavelet analysis on the Chandra data which con- 
firms there is no source detected. Comparison of the 
XMM-Newton field of PSR B0355+54 with the Digitized 
Sky Survey does not show any optical source at the po- 
sition of the south west emission. If the emission in the 
XMM-Newton data is real, it suggests that the diffuse 
emission could be varying over time. 

We smoothed the Chandra ACIS image with a Gaus- 
sian of width ~ 2" (Figure bottom panel), the result- 
ing image suggests that there are two regions of enhanced 
diffuse emission - one near to the pulsar and the other 
~ 10" away. The intensity profiles for the XMM-Newton 
and Chandra data indicate that the core of the X-ray 
emission lies within 5" of the pulsar position. Both emis- 
sion profiles indicate that the diffuse emission extends 
out to ~ 50", with the bulk of the emission lying within 
20 - 30" of PSR B0355+54. The profiles also show ev- 
idence for a dip in the emission at ~ 10" agreeing with 
the smoothed ACIS image. 

Detection of emission from PSR B0355+54 at distances 
of 1.'6 to 5' from the source position i.e. on a lar ger scale 
than shown in Figure HI have been repor ted by Hclfand 
fHWl.lfiewarH Wa.nel ((T988I) andlSktrtel fmol . see also 
iTepedelenlioelu fc Qgelmanl l)2005|) . Visual inspection of 
a smoothed version of the MOS1 image indicates that 
there is a region of enhanced emission extending a few 
arcminutes south of the pulsar, orientated ~ 20° closer to 
North than the diffuse emission we have reported above. 
The intensity of the MOS1 emission in this direction 
was determined using bilinear interpolation at regularly 
spaced intervals from the source. We show in Figure 
the distribution of counts as a function of distance from 
the source compared to the mean background. Our anal- 
ysis suggests that there is an excess of counts at V and 
~ 3 — 5' from the point source. 

4. SPECTRAL ANALYSIS 



In order to investigate the properties of the X-ray emis- 
sion from PSR B0355+54 and the compact (< 50") dif- 
fuse nebula, we compared the spectra extracted from 
different spatial regions. Our results from the spatial 
analysis suggest that the core of the pulsar emission lies 
within 5" of the pulsar's position. However, in the case 
of the XMM-Newton data, this size of aperture does not 
contain enough counts for a meaningful analysis. 

The pulsar spectrum has been extracted from the 
XMM-Newton observation using a circular region of ra- 
dius 30", centered on the pulsar's radio position. The 
background was extracted from a region of similar size 
offset from the pulsar position. The total counts con- 
tained in the source region is 1143 with an estimated 562 
from background. The spectrum was regrouped by re- 
quiring at least 50 counts per spectral bin. We created a 
photon redistribution matrix (RMF) and ancillary region 
file (ARF) for the spectrum. The subsequent spectral fit- 
ting and analysis was performed using XSPEC, version 
11.3.1. 

We modeled the spectrum in the 0.5 — 9.0 keV range. 
Initially we fitted the spectrum with single-component 
models including absorbed power-l aw, blackbody and 
magnetized, pure H atmospheric ijPavlov et alJ Il995(l 
models. The spectrum is best-fitted with a power-law 
with index Y = 1.5^3 and column density Njj = 

(0.50±g;|g) x 10 22 cm" 2 . This value for the power-law 
index is similar to the values found for other PWNe (e.g. 
IKaspi et al.l2005j) . The Galactic hydrogen column in the 
direction of PSR B0355+54 is N H = 0.88 x 10 22 cm" 2 . 
The fit results in an unabsorbed 0.3-10 keV energy flux of 
(2.3±J;?) x 10" 13 ergs cm" 2 s" 1 . We also fitted the spec- 
trum with blackbody plus power-law and atmospheric 
plus power-law models, both modified by photoelectric 
absorption. The multi-component models give similar 
values for reduced % 2 , however the temperatures implied 
by the fits are poorly constrained. It is likely that the 
presence of the pulsar wind nebula, and being unable to 
separate the pulsar core and diffuse emission, effects our 
ability to constrain the thermal component in the spec- 
tral fits of the XMM-Newton data. The results of the 
XMM-Newton spectral fitting are given in Tabled The 
XMM-Newton spectrum with the best-fitting power-law 
model is shown in Figure^ 

For the Chandra data we extracted a spectrum for 
the core of the pulsar emission from a circular region 
of radius 5" centered on the radio position. An annulus 
centered on the pulsar position was used to extract the 
background, with inner and outer radii of 6" and 10", re- 
spectively. We find a total of 244 counts contained in the 
source region, with 29 counts attributed to background. 
We created the RMF and ARF files using standard CIAO 
tools. Before fitting the spectrum we regrouped the data, 
requiring a minimum of 15 counts per spectral bin. 

We fitted the spectrum in the energy range 0.5 — 7.0 
keV using the same models as for the XMM-Newton data. 
In the first instance we let the neutral hydrogen column 
density be a free parameter; however this led to unrea- 
sonably small values for Njj- Subsequently we fixed the 
column density at the value found from the power-law fit 
to the XMM-Newton spectrum of PSR B0355+54. We 
find that the Chandra spectrum can also be character- 
ized by a power-law. The model has a power-law index 
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TABLE 1 

Best-fit parameters for the XMM-Newton spectrum of PSR B0355+54 



Model N H T T/T™ f xl (dof) F x 

(10 22 cm- 2 ) (10 6 K) (erg cm- 2 s" 1 ) 



PL 




U.OU_ 20 






0.7 (16) 


(2.3+J;° 7 ) x 10- 


13 


BB 




< 0.01 




li.26+|;g 


1.0 (16) 


(1.4t?i) x 10_ 


13 


BB + 


PL 


n 9 o+0.69 
U - ZO -0.23 


0.8+° j 


7.66lf 6 1 4 3 


0.7 (14) 


(2.0 + i ?) x 10 ~ 


13 


NSA 




13+ 023 




7.66t££ 


0.9 (16) 


(1.4+°^) x 10- 


13 


NSA 


+ PL 


49+ - 21 


•■•• o -0.3 


0.64+ ;? 3 


0.8 (14) 


(2.4±i;J) x 10- 


13 



Note. — The last column is the unabsorbed flux in the 0.3-10 keV range. In 
the case of the NSA model the mass and radius of the neutron star are fixed at 
Mns = 1-4M0 and -R/vs = 10 km, respectively, and the magnetic field of the 
neutron star is fixed at B = 10 12 G. The errors quoted are the 90% uncertainties. 
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Fig. 4.— The XMM-Newton spectrum of PSR B0355+54 with 
best-fit power-law model. Also shown are the residuals from com- 
parison of the data to the model. 



of r = 1.9+q' 4 . which is consistent within the 90% uncer- 
tainties to the value found from the XMM-Newton data. 
However, in the case of the Chandra data we find that a 
thermal plus power-law model provides a better fit statis- 
tically. The data are equally well-fitted by a blackbody 
plus power-law and a magnetized, pure H atmospheric 
{Pavlov et alTToM "nsa" model in XSPEC) plus power- 
law model. The results of the Chandra spectral fitting 
are given in Table [5] 

For the blackbody plus power-law model the best- 
fit parameters are a power-law index of T = l.O+J'o 
and temperature of T = (2.32±J'|f) x 10 6 K. Us- 
ing a distance to the source of D = 1.04+o' l6 " kpc 
l)Chatteriee et alJl2004(h this implies a blackbody emit- 
ting radius of 0.12±g;Jf km. This value is too small to be 
reconciled with the radius of the neutron star and would 
indicate that the origin of the emission is a hot polar 
cap. We find an unabsorbed 0.3-10 keV energy flux of 
(6.4+. j 9 { 3 ) x 10 -14 ergs cm -2 s -1 . The magnetized, pure 
H atmospheric plus power-law model best-fit parameters 
are a power-law index of T = 1.5+ '4, temperature of 
Tg f = (0.451^2°) x 10 6 K, and a radius for the neutron 
star of Rns = 7.2+2 \ km. For this fit the distance to 



the source and the mass of the neutron star were fixed 
at D = 1.04 kpc and M N $ = 1.4M©, respectively. The 
magnetic field of the neutron star was fixed at B = 10 12 
G (this is a good approximation since the pulsar mag- 
netic field as infer red from radio timing properties is 
B = 8.4 x 10 11 G, iHobbs et"alll2004t Manchester etaD 
120051) . The unabsorbed 0.3-10 keV energy flux for this fit 
is (1.5l^ 4 7 °) x 10~ 13 ergs cm" 2 s" 1 . The Chandra spec- 
trum of the core emission from PSR B0355+54 is shown 
in Figure [3 with the best-fitting blackbody plus power- 
law model (top) and magnetized, pure H atmospheric 
plus power-law model (bottom). 

To analyze the compact diffuse emission we created a 
new events file in which the emission from the pulsar core 
was removed. We extracted a spectrum for the diffuse 
component from a rectangular region of 40" x 55", cen- 
tered on the emission and orientated along the direction 
of the pulsar's proper motion. The background was ex- 
tracted from a region of similar size offset from the diffuse 
emission. The diffuse component extraction region con- 
tains 1207 counts, with an estimated 414 counts due to 
background. We created the RMF and ARF files using 
standard CI AO tools. Before fitting the spectrum we re- 
grouped the data, requiring a minimum of 15 counts per 
spectral bin. We modeled the spectrum over 0.5 — 7.0 
keV with an absorbed power-law, keeping the column 
density fixed at N H = 0.50 x 10 22 cm" 2 . The best-fit 
has a power-law index of T = 1.4 ± 0.3 and unabsorbed 
0.3-10 keV energy flux of (1.7+^) x 10" 13 ergs cm" 2 
s" 1 . 

In order to investigate the possibility of spectral evolu- 
tion along the extended X-ray emission we created spec- 
tra for three regions of the compact diffuse emission. The 
sizes of the regions were chosen with the aim of having 
similar numbers of counts in each region. The three ex- 
traction regions, orientated along the direction of proper 
motion, are as follows, region 1: 40" x 18", contains a 
total of 396 counts with 130 attributed to background, 
region 2: 40" x 9", contains a total of 356 counts with 
69 attributed to background, region 3: 40" x 28", con- 
tains a total of 454 counts with 214 attributed to back- 
ground (see Figure H (middle panel)). The background 
was extracted from the same region as above. We cre- 
ated response files for each region and regrouped the 
spectra, requiring a minimum of 15 counts per spectral 
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TABLE 2 

Best-fit parameters for the Chandra X-ray emission of PSR B0355+54 



Region 




Model 


r 


T/T eff 
/"1 no is\ 
(10 K) 


Rbb/Rns 
km 


XZ (dof) 


Fx 

/„,„ „__ — 2 „-l\ 
(erg cm s J 




Core 




PL 


1 q+0.4 
— 0.3 






0.5 (34) 


(4.9+V 7 ) x 10- 


14 






BB 


6.73_t 2 j 2 


01 + 03 


1.4 (34) 


(2.4«; 1 1 ) x 10- 


14 






BB + PL 


i o +0 - 2 

1.0 


o 09+I.I6 

-0.81 


12+° 16 


0.3 (32) 


(6A+IY) x 10 ~ 


-14 






NSA 




0.48 


9.5 


3.6 (34) 


(l.8±g;|) x io- 


13 






NSA + PL 


i tr+0.5 
L -°-0A 


45+0- 20 


7 9+7.2 
'■ z -2.2 


0.4 (32) 


(Lsj:^ ) x io- 


-13 


Diffuse - 


all 


PL 


1 4+0-3 






1.0 (50) 


(1.7l°'j) x 10- 


13 


Diffuse - 


1 


PL 


1 4+0-4 
0.4 






1.0 (16) 


(5.7l : ^) x 10- 


14 


Diffuse - 


2 


PL 


i =+0.3 
^-OA 






1.4 (16) 


(r.6±?-|) x io- 


14 


Diffuse - 


3 


PL 


i O+0.5 
1 - z -0.4 






1.2 (17) 


(5.3 +4 / 4 ) x 10- 


14 



Note. — The neutral hydrogen column density has been fixed at Njj = 0.50 X 10 22 cm -2 
in all of the fits. The last column is the unabsorbed flux in the 0.3-10 keV range. In the case 
of the NSA model the distance to the source is fixed at D = 1.04 kpc, the mass of the neutron 
star is fixed at Mjvs = 1-4Mq, and the magnetic field of the neutron star is fixed at B = 10 12 
G. The errors quoted are the 90% uncertainties. 
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Fig. 5. — The Chandra spectrum of the core emission of 
PSR B0355+54 with best-fit blackbody plus power-law model (top) 
and magnetized, pure H atmospheric plus power-law model (bot- 
tom). Also shown are the residuals from comparison of the data to 
the model in each case. 



bin. The spectra were fitted in the 0.5 — 7.0 keV en- 
ergy range with a power-law and fixed column density 
of N H = 0.50 x 10 22 cm" 2 . We find the best-fit power- 
law indices for regions 1-3 are T = 1.4 ± 0.4, 1.5 ± 0.3 
and 1.2_tg |, respectively. Due to the uncertainties on the 
indices the presence of any spectral variability remains 
unclear. The unabsorbed 0.3-10 keV energy fluxes are 
(5-7±f:l) x 10- 14 , (7.6±f;|) x 10- 14 and (5.3l£$) x 10~ 14 
ergs cm -2 s _1 , respectively. The Chandra spectra of the 
diffuse emission from regions 1-3 are shown in Figure H3 
(first-third panel) with the best-fitting power-law mod- 
els. 

We have also tried to determine if there are any spec- 
tral changes by using the hardness ratio ft.2.o> which is 
defined as the ratio of counts above 2.0 keV to that be- 
low 2.0 keV. For the whole of the compact diffuse com- 
ponent we find /12 o = 0.97 ± 0.09. Regions 1, 2 and 3 
have h 2 .o = 1-05 ± 0.17, 0.91 ± 0.13 and 0.97 ± 0.18, re- 
spectively. Due to the uncertainties, no particular trends 
in hardness ratio can be determined from one region to 
the next. 

5. TIMING ANALYSIS 

We barycentrically corrected the photon arrival times 
in the XMM-Newton PN event file before performing the 
temporal analysis. We extracted data for the source from 
circular regions of 15" and 30" centered on the pulsar 
position. The total counts encompassed in these regions 
are 391 and 1143 respectively, with the background con- 
tributing 151 and 562 counts, respectively. 

In order to search for an X-ray modulation at the 
PSR B0355+54 spin frequency, we first determined a 
predicted pulse frequency at the epoch of our XMM- 
Newton observations, assuming a lin ear spin-down rate 
and using the radio m easurements ( Ho bbs et al.1 12001 
Manchester etHl 12005 ) . We calculate / = 6.3945388 
Hz at the midpoint of our observation (MJD 52315.7). 
As glitches and/or deviations from a linear spin-down 
may alter the period evolution, we then searched for a 
pulsed signal over a wider frequency range centered on 
/ = 6.39454 Hz. We searched for pulsed emission using 
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Fig. 6.— Spectral fitting to Chandra data of PSR B0355+54. 
First-third panels: diffuse emission from regions 1-3, respectively, 
with best fit power-law models. Also shown are the residuals from 
comparison of the data to the model in each case. 



two methods. In the first method we implement the Z% 
test (jBuccheri et al.lll983(l , with the number of harmon- 
ics n being varied from 1 to 5. I n the second method we 
calcu late the Rayleigh statistic (Ide Jaeerl 119911 iMardial 
Il972j) and then calculate the maximum lik elihood peri- 
odogram (MLP; see e.g. IZane et all 12002^1 using the C 
statistic i|Cashlll979]) to determine significant periodici- 
ties in the data sets. 
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Fig. 7. — Z 2 -test (top) and Maximum Likelihood Periodogram 
(MLP; bottom) for the PN data of PSR B0355+54. The dom- 
inant peak from the Zj-test and the corresponding peak in the 



MLP are marked. The peaks occur at 6.3945447 



.3945467 



+0.00000821 



+0.0000167 
-0.0000107 



Hz and 



-0 00000818 res P ec ti ve ly- The dotted lines represent 
the 68% (x 2 = 2.71) and 90% (x 2 = 1.0) confidence levels for the 
frequencies in the MLP. 



The frequency search of the data extracted from an 
aperture of radius 30" does not yield any significant 
peaks near to the predicted frequency with either search 
method. From our spatial analysis we know that the core 
of PSR B0355+54's emission lies within ~ 5" of the pul- 
sar position. We have therefore also searched for pulsed 
modulations in data extracted from a smaller aperture, 
however we have used a radius of 15" as any smaller does 
not encompass enough counts for a meaningful analysis. 

The most significant Z^-statistic occurs for n = 1. 
With the number of harmonics equal to one, the Z\- 
statistic corresponds to the well known Rayleigh statis- 
tic. We find three peaks with > 90% significance in the 
MLP, all with corresponding peaks from the Zf-test (see 
Figure [TJ. The dominant peak from the Zf-test occurs 



at 6.3945447 



+0.0000167 
-0.0000107 



Hz, with the corresponding peak 
in the MLP occurring at 6.3945467±fJ .qoooosis Hz - The 
uncertainties quoted are the 68% confidence limits on 
the position of the peak. Both frequencies are consis- 
tent, within the 68% contour, with the predicted pulse 
frequency, and with each other within the 90% contour. 
The second most prominent peak from the Z^-test, and 
the corresponding peak in the MLP, are not consistent 
with the predicted pulse frequency. 

While wc have detected a frequency that is in 
agreement with the predicted pulse frequency for 
PSR B0355+54 we caution that the Z\ peak has a prob- 
ability of chance occurrence of 3 x 10~ 3 . Further observa- 
tions of the source are needed to show whether the mod- 
ulation detected is in fact pulsed X-ray emission from 
PSR B0355+54. We have folded the data on the pre- 
dicted pulse frequency and the frequency found from the 
Z 2 -test (see Figure EJ; by fitting the profiles with a si- 
nusoid wc find that the modulation amplitude for the 
former is 25 ± 7% and 21 ± 8% for the latter. 

6. DISCUSSION 

Our spatial analysis of the XMM-Newton and Chandra 
observations of PSR B0355+54 have not only revealed X- 
rays from the pulsar, but have provided definitive proof of 
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Fig. 8. — PN data in the 0.3-10 keV energy range for 
PSR B0355+54 folded on the frequency predicted from the radio 
measurements (top) and the frequency found from the Z^-test (bot- 
tom). In both cases the data are folded using the radio ephcmeris. 



diffuse emission extending in the opposite direction to the 
pulsar's proper motion. Similar detections of extended 
emission have been see n for other sou rces (e.g. N157B, 
Wang fc GottIell[l998t PSR B1757-24. lF7ail fc Kulkarrii 
1991 iKasni et alMOOlt PSR B1957+2 0, lStappers et al. 
20031 PSR B1951+32. [Li et al]l2005l) . and have been 
interpreted as emission from a ram-pressure confined 
PWN. 

We cannot separate the core and diffuse emission com- 
ponents for the XMM-Newton data and find that the 
spectrum can be well-fitted with a power-law model with 
index T = 1.52.°;^, simila r to the value found for other 
PWNe i|Kaspi et al.l2005r) . The nebular emission is most 
likely dominating the spectrum. The core emission from 
the Chandra data can be well-fitted by a thermal plus 
power-law model. A fit with a blackbody plus power-law 
model gives T = (2.32± ;^) x 10 6 K and T = 1.0±°%. 
The fitted blackbody flux corresponds to an emitting ra- 
dius of O.^Iq'Jy km, where the distance to the source 
is 1.04 kpc. In this case, the size of the emitting re- 
gion implies that the flux originates from a hot polar 
cap. We can also fit the spectrum with a pure H, mag- 
netized atmospheric plus power-law model. With the 
distance to the source fixed at 1.04 kpc, this fit re- 
sults in T£ff = (0.45t{j;^) x 10 6 K, T = 1.5±g;| and 



Rns = 7.2^2 2 km. Taking into account the possible 
detection of pulsed emission from PSR B0355+54 it is 
likely that the emitting region is a hot polar cap. 

It is suggested that the presence of a PWN is re- 
lated to the spin-down power of the pulsar, and for 
sources with \ogE < 3 6 the PWN emission efficiency 
is significantly reduced (iFrail fe Scharringhausen|ll997t 
iGaensler et all200(llOotthell200fl . For PSR B0355+54 
log-E = 34.6, making it one of a handful of sources with 
spin-down pow er below this limit wit h a detectable PWN 
(cf. Geminga, iCaraveo et alJl2003T) . Using the results 
from the blackbody plus power-law fit to the core emis- 
sion detected with Chandra we determine an isotropic un- 
absorbed luminosity in the 0.3-10 keV band of 8.3 x 10 30 
erg s _1 . With E = 4.5 x 10 34 erg s" 1 for PSR B0355+54, 
this leads to a conversion efficiency of 2 x 10~ 4 . So in 
fact, we find that the conversion efficiency of the point 
source is similar to the values fo und for other pulsar s 
(see e.g. iBecker fc Trumped IT9971 IGaensler et al.ll2004j) . 
Our analysis also indicates that the compact diffuse com- 
ponent is more luminous than the point source, with 
a conversion efficiency of 5 x 10~ 4 in the 0.3-10 keV 
range. This result is again consistent with other sources 
l)Becker fc Trumped I1997T) . In addition, it is reported 
that when the pulsar spin-down energy is \ogE < 36.5 
the morphology of th e PWN seems to transition from 
toroidal to a jet/tail ijKaspi et al.ll2005|) . Our measure- 
ments of PSR B0355+54 appear to agree with this trend. 

The morphology of the diffuse emission depends on 
how the interaction with the interstellar medium (ISM) 
or s upernova remnant constrain s the flow of particles 
(e.g. [Reynolds fc Ohevalierlll984[) . For a pulsar that is 
moving with a supersonic space velocity, the interaction 
of the supersonic flow with the ambient medium causes 
the speed of the flow to decrease sharply, while the den- 
sity increases, forming a bow shock. In addition to the 
bow shock, which is at some distance ahead of the pul- 
sar, a reverse shock is formed nearer to the source which 
terminates the pulsar wind. 

The results of the spatial analysis of the XMM-Newton 
PN and Chandra ACIS data of PSR B0355+54 indi- 
cate that the bulk of the diffuse emission extends ~ 50" 
[0.25(g?/1.04 kpc) pc] downstream from the pulsar. Us- 
ing the measurements of the pulsar's proper motion 
l)Ghatteriee et all 120041) we find that the transverse ve- 
locity of PSR B0355+54 is v t = 61 km s _1 . This implies 
that the time taken for the pulsar to have traversed the 
length of the diffuse emission is > 4000 yr. In addi- 
tion, by considering the analysis of the XMM-Newton 
MOS1 data we find that the diffuse emission may ex- 
tend as far as ~ 5' [1. 51 (d/ 1.04 kpc) pc] from the point 
source. This results in a travel ti me of > 24000 yr 
for th e pulsar. Following the wor k of iWang fc GotthelJ 
(1998) (see also iKaspi et al.ll2001|) the synchrotron life- 
time of an electron of energy E (in keV) can be defined 

as t s ~ A0E~ l l 2 B_ S J 2 yr, where B_^' 2 is the magnetic 
field in units of 10 -4 G. Assuming that the dominant 
loss mechanism is synchrotron emission, i.e. B > 3.2 
fj,G, and that the energy of the photon is E ~ 5 keV, 
then t s ~ 3000 yr. This indicates that the diffuse emis- 
sion that we detect is not due to particles deposited by 
the pulsar as it traveled through space. Hence, there 
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must be a constant supply of wind particles traveling at 
velocities greater than the space velocity of the pulsar. In 
addition, the particle flow velocity must be high enough 
such that the time for the flow to cross the length of the 
diffuse emission is less than the radiative lifetime of the 
particles. 

Using the Chandra data we have modeled the spectrum 
of the compact diffuse emission, excluding the contribu- 
tion from the pulsar, finding that the data can be well- 
fitted with a power-law. In other sources the power-law 
is seen to so ften as one moves away from the pulsar posi- 
tion (see e.g. lSlane et al.l2002tlLi et all2005l:lKasT)i et all 
2005). An increase in the spectral index is expected as 
the particles will be cooler, i.e. older, at greater dis- 
tance from the pulsar. Our results indicate that we are 
detecting relatively hard emission, but due to the uncer- 
tainties, we are unable to comment on any changes in the 
spectral slope. To measure cooling the PWN must be of 
an adequate size, it may be that for PSR B0355+54 the 
compact diffuse region is not large enough for a substan- 
tial change in power-law index to be measured, and there 
are too few counts in the more extended diffuse region to 
perform a spectral analysis. It is noted however that by 
comparing the spectral indices from the blackbody plus 
power-law fit to the core emission and the power-law fit 
to the compact diffuse emission we do detect an increase 
i n r of ~ 0.5. 

IGaensler et al.l l)2004|) have presented a detailed analy- 
sis of the diffuse X-ray emission associated with the radio 
source G359.23-0.82, also known as "the Mouse". Their 
hydrodynamic simulations show that there are a num- 
ber of regions that can be defined in a pulsar bow shock. 
These include a pulsar wind cavity, shocked pulsar wind 
material, contact discontinuity (CD) and shocked ISM. 

The energetic shocked particles from the pulsar are 
confined by the CD, the position of which denotes the 
transition to the sho cked ISM. Following the method of 
IGaensler et"afl <|2004) we have estimated the distance be- 
tween the peak of the emission from PSR B0355+54 and 
the sharp cut-off in br ightness ahead of the pulsar. Us- 
ing the same limit as IGaensler et alJ l)2004l) i.e. where 
the X-ray surface brightness falls by 1/e 2 = 0.14, we 
find a distance of 0.9" ± 0.2", giving the CD a projected 
radius of tcd = 0.004 ± 0.001 pc. Here, and in the fol- 
lowing, we have used a distance to the pu lsar of 1.04 kpc 
llChatteriee et al.ll2004|) . From Eq. (1) of IGaensler et all 
( 2004) we can estimate the radius of the forward termi- 
nation shock (TS), r^ s ~ 0.003 pc. This corresponds to 
an angular distance of 9 = 0.59". Comparing our val- 
ues to those for the Mouse implies that the emission in 
front of the pulsar is more compact in PSR B0355+54 
than for the Mouse. In both cases the close proximity 
of the forward TS to the peak X-ray emission renders 
the TS undetectable. Using our results and Eq. (2) of 
IGaensler et all p00|) we find that PSR B0355+54 pro- 
duces a ram pressure of pv\ ~ 1.4 x 10 -9 ergs cm -3 . 

— 1/2 

Assuming cosmic abundances, this gives Vt <~ 247n 
km s , where tiq is the number density of the ambient 
medium, so for PSR B0355+54 we determine no ks 0.06 
cm~ 3 , which is not unrealistic. 
Additional information can be obtained by equating 



the pressure of the pulsar wind (assumed isotropic), to 
that of the ambient medium. By introducing the Mach 
number M = v t /c s , where c s is the adiabatic sound speed 
in the ambient medium, and using the same prescription 
as IGaensler et alJ l)2004D for a representative ISM pres- 
sure (i.e. Pism = 2400fcP erg cm -3 , with 0.5 < P < 5 
and k is the Boltzmann's constant), this gives: 

E/iMrrsfc] = 2400fc 7/S A/iW 2 , (1) 

from which we can obtain an estimate of the Mach num- 
ber. We find that for PSR B0355+54 the sound speed c s 
of the medium lies in the range 1 — 30 km s . The three 
principal phases of the ISM are generally named cold, 
warm and hot and are characterized by typical sound 
speed values of 1, 10 or 100 km s _1 ; according to this de- 
nomination our result implies that the pulsar is moving 
in either a cold or mildly warm a mbient gas. For com- 
parison, in the case of the Mouse, Gacns ler et al.l l)2004fl 
found that the most probable pulsar velocity requires 
that the pulsar is moving through a warm phase of the 
ISM 

IGaensler et alJ l)2004f) also discuss the possible detec- 
tion of the backward TS in their data. Their simulations 
show that this feature has a closed surface, while the 
CD and bow shock are unrestricted. The backward TS 
should lie much further away from the pulsar than the 
forward TS, i.e. r^ s 3> In principle this means that 
the backward TS may be detectable. The possible dip 
we see in the profiles for the PSR B0355+54 data could 
indicate the presence of the backward TS. The angular 
separation of the dip in our data is ~ 10", a value con- 
sistent with that for the Mouse. For a backward TS, 
IGaensler et alJ l)2004|) predict that there would be a lack 
of spectral evolution, a result we have found for the dif- 
fuse emission of PSR B0355+54. However, we note that 
the feature in the Mouse data (and simulations) is quite 
compact in the north-south direction in comparison to 
the PSR B0355+54 feature. In addition, the number 
of counts we detect for PSR B0355+54 may hinder our 
investigation of the presence of such a feature. Deeper 
observations are needed to probe further the PWN of 
PSR B0355+54. 

A few days before submission of this paper, Tepedelen- 
lioglu & Ogelman submitted a paper based on the same 
(public) observations to ApJL (astro-ph/05 12209). 
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